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The  problem  of  standing  and  progressing  waves 
on  a  beach  with  a  plane  bottom  has  been  solved  by 
Miche  [1],  Lewy  and  Stoker  [2],    [3]  for  slope 
angles  a>  which  are  integral  fractions  of  a  right 
angle,'""  co  «=  rr/Sn,  n  =  1,2,,..  •  In  the  present 
report,  in  which  various  suggest ions  by  J«  J,  Stoker 
and  H.  Lewy  are  incorporated,  a  simple  asymptotic 
representation  for  small  angles  co  la  presented. 
It  enables  us  in  particular  to  show  that  for  small 
slope  angles  the  flow  can  be  described  as  consist- 
ing of  waves  with  a  varying  wave  length,  which 
depends  on  the  depth  in  the  same  way  as  Airy's 
waves  in  channels  with  constant  depth.  In  addition 
our  formulas  give  the  dependence  of  the  amplitude 
and  the  phase  on  the  depth  for  small  slope  angles. 
The  wave  shape  calculated  on  the  basis  of  our  formulas 
for  a  beach  with  a  6°  slope,  see  Figures  2  and  3, 
agrees  almost  perfectly  with  that  calculated  on  the 
basis  of  the  exact  theory,  see  Stoker  [3,  p,  35  ] , 


*Lewy  solves  more  generally  the  problem  with  the 
slope  angle  <^    =  pTr/2n,  p  =  1,3,5,,.,   ,  (p,n)  =  1. 


The  derivation  of  this  asymptotic  representation 

from  a  new  form  of  the  solution  for  slope  angles 

«^  =  Tr/2n  will  be  given  in  a  later  publication. 
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1,  Mathematically  speaking  the  problem'"  is 
that  of  finding  an  analytic 
fijnction   X(z)  of  the  com- 
plex variable  z  =  x  +  ly 
defined  in  the  sector 

-  ^  <  ©  <  o» 


z   =  rel*. 


r  >  0.  see  Figure  1  .  which 


_.      * 


satisfies   the  bovmdary  con- 
ditions 


VJoler    turfac* 


•T^ 


bottom 


pig-ore  1 
Sloping  beach 


In  i-S-   +  1)  X  =0   on  the  line  ©  «  0, 


II 


Im  e"^^  ^  ^  ~  °   on  the  line  9  =  -  o;  , 


remains  boimded  as  r  approaches  infinity,  and  be- 
comes at  most  logarithmically  infinite  as  r  -*•  0, 
The  function 


(1) 


)C(z,t)  =   %(z)cos  t 


is  the  complex  potential  function  of  the  flow  rep- 
resenting a  standing  wave.  The  units  of  z  and  t 
are  so  chosen  that  both  the  circular  frequency  5" 
and  the  acceleration  of  gravity  g  have  the  value 
one,  see  Stoker  [3,  Section  3],  The  velocity 
coinponents  u  and  v  in  the  x  and  y  directions  are 


The  problem  of  waves  on  sloping  beaches  and  the 
connection  of  the  questions  here  treated  with 
related  questions  are  explained  in  an  introductory 
report  [4], 
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derived  from  the  function  X  through 
u  -  Iv  =  dX^dz.  The  boundary  condition  II  ex- 
presses the  requirement  that  the  water  remain 
attached  to  the  rigid  bottom  surface  0  =  -  t*^  • 
Condition  I  resxxlts  from  linearizing  the  condition 
that  the  presstire  on  the  top  surface  remain  con- 
stant. We  note  that  the  elevation  of  the  free 
water  surface  is 

(2)         t^U,t)  =  Re  ^   X(x,t). 

For  details  see  the  two  papers  [2  and  3]  quoted, 

2,  There  are  two  solutions  of  the  problem. 
One  solution,  X  (z),  is  regular  at  the  origin, 
z  -  0,  and  behaves  like   X?,  =  -iC  e"^^  at  infinity 


00 


(3)    X^(z)  +  IC  e"^^  -s-  0   as   |z|  ->  oo 


with  an  appropriate  complex  constant  C,  see  [3,  Sec.  2], 
The  other  solutionj 
in  such  a  way  that 


The  other  solution,   X^^z),  is  singular  at  z  =  0, 


(4)     X^(^)  "  2  ^°S  ^   ^^  regular  at   z  =  0, 


while   X^  behaves  like  X^  =  C  e"^^  at  infinity. 


(5)     X^(z)  -  C  e"^^  ->  0       as  |z|  -*  00 


-4- 

The  standing  waves,  given  by 

yR  =   X^(z)cos  t   and     X^  =   X^(z)sin  t  , 

ape  evidently  90°  out  of  phase  at  infinity. 

With  the  aid  of  these  two  standing  waves  one 
can  build  up  progressing  waves  traveling  toward 
shore.  Such  a  progressing  wave  is  given  by  the 
complex  potential 


(6)     X^(z,t)  =  X^(z)cos  t  +   X^(z)sin  t  , 


which,  at  infinity,  behaves  like 

(7)     X^(z,t)  -  c  e-'^^-*-^^  z  ;soo  , 

and  thus  like  a  sinusoidal  wave  traveling  toward 
shore . 

In  the  following  it  is  convenient  to  introduce 
the  functions 


X"(2)  =^(X^(z)  +  iX^(z)). 


V/aves  progressing  toward  shore  are  then  given  by 


(9)  X^(z.t)    =      lCiz)e^^  +      X"(z)e-"      . 


3,   Ovir  main  problem  is  to  find  an  asymptotic 

representation  of  the  functions    X  (z)  and 
g 

X  (z)  for  sinall  values  of  the  slope  angle  co  • 

There  are  tliree  different  asymptotic  representations 

depending  on  which  quantities  are  kept  fixed  while 

one  lets  co  approach  zero. 

The  first,  aiid  raost  natural,  problem  of 

asymptotic  representation  arises  when  one  asks  for 



the  asymptotic  value  of   X  (z)  aJ^d    X.  (z)  at  a 

fixed  point  z.  Evidently,  only  points  z  =  x  on  the 
water  surface  remain  in  the  water  region  for  arbi- 
trarily small  angles  i*^  ;  therefore,  we  consider 
only  points  z  =  x.  In  accordance  with  the  shallow 
water  theory,  see  Stoker  [3,  Section  3],  the 
answer"  to  this  question  is  probably  that  X  (x) 
and   X   (x)  behave  asymptotically  like 


(10)  Xi  =  2TrJ^(2  \/v^  ),  11  =   SttYqCS  \f^7^  ) ,      co   Z  0   , 


involving  the  regular  and  singular  Bessel  f\anctions 

J^  and  Y  •  These  Bessel  functions  die  out  at  infinity, 
o      o 

Therefore  the  approximate  solutions  (10)  cannot 
represent  the  exact  solutions  far  away  from  shore 
since  the  wave  motion  does  not  die  out  but  beloaves 
sinusoidally  at  large  distances  from  shore.  For  a 
slope  angle  <^  =  6   ,   for  example,  the  expressions 


'The  approximate  solution  for  beaches  with  a  small 
slope  based  on  a  power  series  development  of  the 
exact  solutions,  given  by  Miche  [1,  p.  139],  is 
essentially  the  same  as  that  given  by  the  shallow 
water  theory. 
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(10)  give  a  good  approximation  only  in  a  region 
covering  one  to  two  wave  lengths  near  shore,  while 
the  asymptotic  representation  to  be  discussed  next 
gives  nearly  perfect  agreement  with  the  exact 
solution  all  along  the  water  surface  except  on  a 
section  near  shore  covered  by  about  one  third  of 
the  length  of  the  first  wave,  see  Figures  2  and  3. 

4,   In  the  second  asymptotic  representation, 
which  we  shall  derive  in  great  detail,  the  product 
O)  z  is  kept  fixed  while  co  approaches  zero.   In 
other  words,  we  ask  for  the  asymptotic  values  of  the 
functions   X  (z)  at  places  z  which  move  away  from 
shore  so  that  their  distance  from  shore  is  inversely 
proportional  to  the  angle  CO    ,  Again  we  consider 
only  points  z  =  x  on  the  water  surface.  The 
asymptotic  representations   X  t   ^^^   "^2  ^-^   ^ 
and   X  "  are  given  by  the  expressions 


+  +iCO"^k(A)+i-? 

(11)    X:(x)  =  -(1  +  X)   A(A)e-        -  ^ 


in  which  A  is  a  function  of  to  x  given  by  the 
inverse  of  the  relation 

(12)   CO  X  =   AArtanhA  =^    Aln  -ji^  ,   0  <,  A  <  1. 
A(X)  is  given  by 

and  k(A),  j(A)  are  given  by 
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(14)         k(A)  =  ArtanhA  +  /  ArtanhV  ^ 

0  ^ 


(15)         j(X)  =  (l-A^)>k»(A)  =  A  +  (1- A.^)Artanh  A 


The  beha,vior  of  the  quantities  k.  A,  J,  and  A/  \/  Sir  to 

as  functions  of  co  x  is  shown  in  Figures  4,  5,  6, 

and  ?•  The  asymptotic  representations  of  the 

R  S 

complex  potential  functions   X  (x)  and   X  (x)  for 

standing  waves  are  then  derived  from  (8),  The 

asymptotic  representation  of  the  complex  potential 

function  %    {x,t)   for  progress: 

toward  shore  is  found  from  (9), 


p 
function  %    {x,t)   for  progressing  waves  traveling 


^  l[co-lk(A)+t+5] 

(16)  tlU,t)    =   -  (1-  A)A(A)e 

-l[w"^k(A)+t+5] 
-  (1+  A)A(A)e  * 


Each  of  the  two  contributions  here  evidently 
represents  a  wave  traveling  toward  shore  without 
being  a  solution  itself. 

The  behavior  of  the  asymptotic  solutions  (11) 
at  infinity^  CO  x  s  oo ,  is  determined  from  the 
formulas 


(17)  X^l  -  2e'^^^   , 


(18)  k(X)  ?:  -^  +  coTi,      j(Aj  »  1  , 
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(19)  A(A)~  \j2^tUi      , 

which  can  be  derived  from  (12),  (14),  (13),  and  (15), 
Hence  we  have  for   co  x  si  oo 

(20)  Xj(x)  2J  0,   Xr(x)  7i    2C  e"^   , 


with 


(21)        C  =  -  \J2^XVM    e 


r2 


-i(co-ll^  +^) 


whence 

(22)  X\{yi)-    X^=  -iC  e-^,  X|(x)»  X^=  C  e"^  , 

as  CO  X  ss  00  , 

while  for  the  progressing  wave  given  by  (16), 

(23)  X2(X't)~  X^  =2C  e"^^''"^*^  as  u>x«0Q. 

The  behavior  of  the  asymptotic  solutions  (11) 
for  small  values  of  Cu  x  is  determined  from  the 
formulas 

(24)  A  «     \fa75     , 

(25)  k(A)   ~     2   \fSix  ,      j(A)J5J    2   \fZox  , 
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(26)  A(7l)  2i     VW    ^w/x      , 

which  are   also  easily  derived  from   (12),    (14),    (15), 
and   (15),     Hence  we  have  for       co  x  ^  0, 

+  . +21  n/x/C     +  iTr/4 

(27)  Xi(x)  2J    -  V^    V  ""/^  Q"  "  ' 


(28)  X|(x)-    2  VtT    '^J'u/Sl  3ln{2Vx/ui    +  Tr/4) 
■K|(x)«    -2  yif    "^fW^  cos(2Vx/w    +  Tr/4), 

(29)  X2(x»t)-     -2  Vr  ^^/T57x  cos(2VxA>     +  t   +  Tr/4) 


T5 

We  note  that  both  asymptotic  solutions   Xo  ^^'^ 
Xo  have  at  to  x  =  0  the  same  singularity,  that  of 


#^, 


which  thus  overshadows  the  logarithmdLc 


singularity  of  the  exact  solution  X     and  the 
regularity  of  the  exact  solution   X  •  The  reason 
for  this  discrepancy  is  that  the  accuracy  with  which 
the  asymptotic  solutions  (11)  represent  the  exact 
solutions  decreases  as  x  approaches  zero  for  a  fixed 
value  of  uj  while  it  increases  as  co     approaches  zero 
for  a  fixed  value  of  co  x. 

It  is  quite  llliominating  in  this  respect  that 
the  functions  (28),  representing  the  "second"  asymptotic 
solutions  (11)  for  small  values  of   cu  x,  at  the  same 
time  represent  the  "first"  asymptotic  solutions  (10) 
for  large  values  of  x,  as  follows  from  the  well  known 
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asyn^Jtotic  expansions  of  the  Bessel  f\inctions  for 
large  values  of  the  argument.  For  a  slope  angle 
^*^  =  6°  the  expressions  (28)  represent  the  solutions 
quite  well  in  a  range  covering  two-thirds  of  the  first 
wave  and  most  of  the  second  wave,  see  Pigiores  2  and  3, 
In  this  range,  therefore,  the  "first"  and  the  "second" 
asymptotic  representations  agree. 

The  s\3rface  elevations  of  the  standing  waves 
are  by  (1),  (2),  (8),  and  (11)  given  asymptotically 
as 

(30)  ^^(x,t)  =  -2A(X)sin(co-lk(»  +  |)sint 

v\|(x,t)  =  2A(A)cos(Cu»"-'-k(A)  +  |)sint. 

As  "local"  wave  length  for  these  surface  waves 
we  may  introduce  the  quantity 

(31)  A=  2Tr[4(co-lk(A)  +  |)]-1  =  2Tr  ^/^ 

because  this  expression  would  be  the  wave  length  If 
^''  k(A)  were  proportional  to  x  and  the  wave  accordingly 
purely  sinusoidal.  Prom  (12)  and  (14)  we  find 
doJx/dA  =  ^k»(A);  hence  the  local  wave  length  is 

(32)  A  =  2TrA  . 

The  coefficient  2A(X)  of  the  sinusoidal  functions  in 
the  expressions  (30)  for  the  surface  elevation  may 
be  called  the  "amplitude"  as  fimction  of  the  wave 
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length.   It  is  remarkable  that  this  an^litude,  which 
according  to  (19)  approaches  the  value  /'StToJ  as 
6Jx  ->•  00 ,  decreases  as  oOx   decreases  until  it 
reaches  the  minimum 

(33)  2A^n  =  '^^^    \^   ^*   cJx  =  1, 

and  then  increases  indefinitely  as  CO  x.   approaches 
zero,  see  Figure  7»  This  fact,  already  observed  by 
Rankine,  see  Micho  [1,  p,287]  is  easily  derived  from 
(12),  (13),  and  (15).   It  is  in  agreement  with  the 
behavior  of  the  exact  solution,  see  [3]» 

For  the  surface  elevation  of  the  progressive  wave 
we  obtain  from  (16)  and  (2) 

(34)  -filiXft)   «  2A(X)sin(a;-lk(\)  +  t  +  |) . 

As  for  standing  waves  we  can  interpret  the  quantity 
2Tr  A   as  the  "local"  wave  length  and  2A(A)  as  the 
"amplitude"  as  function  of  the  wave  length.  The 
expression  (13)  for  this  amplitude  agrees  with  that 
derived  by  Miche  p.,p2851  for  progressing  waves  on  the 
basis  of  energy  flxix  considerations, 

5,  We  now  proceed  to  discuss  the  third  asymptotic 

problem.   It  arises  if  one  considers  as  fixed  the 

distances  in  x-direction  from  that  cross  section  x  =  x^ 

at  which  the  depth  of  the  water  has  a  prescribed  value 

h  •  Thus  the  shore  line  recedes  to  -oo  as  to  — >  0, 
o 

If  the  angle  cJ    is  small  we  may  set  a)x  =  h.  Hence 
the  cross  section  x  =  x  is  given  by 
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(35)  Xq  =  ^"\  . 

Accordingly,  we  keep  the  values  z  -  x  =  z  -  co*''Ti 

fixed.  The  asymptotic  expressions  of  the  third  kind 

for  the  functions  \  iz)   at  the  top  z  =  x  of  the 

cross  section  x  =  x  are  obviously  just  the  saine  as 

the  asymptotic  expressions  of  the  second  kind.  The 

asymptotic  expressions  of  the  third  kind  at  a  point 

z  =  X  with  fixed  x  -  x  can  also  be  derived  from  those 

of  the  second  kind;   for,  a  fixed  value  of  x  -  x 

corresponds  to  values  of  co  x  -  c«j  x  =  CtJ  x  -  h^ 

o  o 

which  approach  zero  as  a)  approaches  zero.  Hence 
the  first  terms  of  the  expansion  of  the  expressions 
of  the  second  kind  with  respect  to  powers  of 
Wx  -  h  give  the  desired  asymptotic  expressions 
of  the  third  kind.   By  virtue  of 


(56)  ^  ^(^)  =  ^"^  ' 


see  (31)  and  (32),  we  have 


(37)      k(X)  =  k(  Aq)  +  Aq^(u;x  -  h^)  +  ...  , 


the  value   A_^  being  determined  from 


(38)  h^  =   A  ArtanhA^ 

o      o       o 


according  to  (12)  and  (35),  Consequently,  we  find 
from  (11)  the  expression 


+lco-^k(A^)+i5  +iX-l(x-x^) 


(39)  X:(x)  =  -(1  +  A^)A(  \)e'  o  -  4  .-  o 


e 
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as  asymptotic  representation  of  the  third  kind  of 
the  functions   X.  (z)  at  the  top  surface  z  =  x.   It 
is  plausible  that  the  asymptotic  representation  of 
the  third  kind  all  through  the  water  depth  near 
z  =  X     is  given  by 

+  +iw-^k(A„)+i5  H-i^l'-'-Cz-x^) 

(40)  XgCz)  =  -(1  +  \)A(Ao)e-       o  -  4  ^-  o     o 


-^o  -  y  ^  °  • 


For  the  standing  waves  we  find 


(41)  7e|'^(2)  =  (l-A.)C^'^  e  °  ^"^° 


-iA:^(2-x  ), 


+  (1+ A^)C^'^  e    °     ° 


valid  for  -h  <  y  <  0,  with 


(42)  c^  =  Cq  ,   C^  =  iC^  , 


ico-lk(  A^)+i? 


(43)  C^=A(X^)e         °   * 


The  standing  waves  (41)  are  evidently  superpositions 
of  two  standing  waves  with  the  wave  length  2Tr  A  l  as 
a  matter  of  fact,  they  are  just  Airy«s  solutions  for 
the  flow  in  a  channel  of  constant  depth  h  .  They 
satisfy  the  conditions 
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(44)      im  (4  +  1)  Xg  =  0 

d 


on 


^^  ^  ^3  =  °     °^     y  =  -^o  ' 


as  is  easily  verified.  Relation  (38)  is  just  AiiT's 
relation  between  wave  length.  Str  A^  and  depth  h  • 

The  question  has  frequently  arisen  whether  or 
not  it  is  correct  to  employ  Airy^s  relation, derived 
for  various  constant  depths^  to  describe  also  the 
variation  of  the  wave  length  on  a  sloping  beach, 
Otir  result  shows  that  it  is  indeed  correct  to  do  so 
approximately,  but  with  good  accuracy,  for  small  slope 
angles. 

The  asymptotic  expression  of  the  third  kind  for 
progressing  waves  traveling  toward  shore  is  derived 
from  (8),  (9),  and  (41) 

P, it^o^(^-^o)^^] 


(45)    XJt(z,t)  =  -(1  -  -KjG^e 


3 


_  -i{x-l(z-xj+t] 


+  (1  +  Ao)C^e 


Both  contributions  are  again  waves  traveling  toward 
shore,  thoiigh  not  solutions  individually.  The  differ- 
ence in  the  two  contributions  lies  in  this  -  that  the 
amplitude  of  the  first  contribution  X 3(2)0   increases 
like  e"^  from  the  top  surface  to  the  bottom  while  that 
of  the  second  contribution  X 3(2)0"   decreases  like 
e^.  Note  that  for  infinite  x  the  first  contribution 
disappears  and  the  second  one  occiors  with  the  amplitude 


-15- 


2  VSttcj,  see  (20),  (21),  while  at  the  shore  both 
contributions  occvlp  with  equal  amplitude,  see  (27). 
The  expression  (45)  of  the  progressing  wave  with 
|C  I  =  A(  Xq)*  except  for  the  dependence  of  the  phase 
on  the  depth,  contained  in  the  expression  (43)  for 
C  ,  was  found  by  Miche  [  1  p,28S]  on  the  basis  of 
considering  the  energy  flux  through  any  cross  section. 
Our  investigation  now  shows  that  this  formula  is  the 
correct  asymptotic  expression  for  small  values  of  the 
beach  angle. 
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